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Abstract. We investigate the occurrence of Hund's rule magnetism in Cg(f molecular ions, by computing 
the ground-state spin for all charge states n from —3 to +5. The two competing interactions, electron- 
vibration (e-v, including Jahn Teller, favoring low spin) and electron-electron (e-e, including Hund-rule 
exchange, favoring high spin), are accounted for based on previously computed ab-initio coupling parame- 
ters. Treating the ion coordinates as classical, we first calculate and classify the static Jahn- Teller distorted 
states for all n, inclusive of both e-v and e-e effects. We then correct the adiabatic result by including the 
zero-point energy lowering associated with softening of vibrations at the adiabatic Jahn- Teller minima. 
Our overall result is that while, like in previous investigations, low-spin states prevail in negative ions, 
Hund's rule high spin dominates all positive Cg^ ions. This suggests also that Hund-rule magnetism could 
arise in fuUerene cation-based solid state compounds, particularly those involving CgJ. 

PACS. 36.40.Cg Electronic and magnetic properties of clusters - 61.48.-l-c FuUerenes and fuUerene-related 
materials (structure) - 71.20.Tx FuUerenes and related materials; intercalation compounds (electronic 
structure) - 75.75.-|-a Magnetic properties of nanostructures 



1 Introduction 



Magnetism without transition metals is a potentially ex- 
citing subject. As an example, carbon magnetism has gen- 
erated some recent interest, in connection with some ful- 
lerene-derived carbon-only magnetic materials |2E|- In 
these materials however the fullerene cage is disrupted into 
some sort of three-dimensional bonding network. In this 
paper we focus on carbon magnetism that may occur in 
isolated fullerene ions, and in ionic fullerene compounds 
where electrons are added or subtracted to Cgo molecules 
that preserve their overall molecular integrity. 

Unconventional properties, including magnetism, of neg- 
ative Geo ions have in fact been previously discussed in 
the literature. Due to the high symmetry of the molecule, 
the threefold degenerate tiu molecular orbital of Geo will, 
when partly filled, be affected by Goulomb exchange |30 
H], which favors molecular Hund's rule magnetism, quite 
similar to that leading to the atomic magnetic moments in 
ordinary d and / elements and compounds. In addition to 
displaying Hund-rule physics, a high-symmetry molecule 
will however also undergo Jahn- Teller (JT) distortions. 
The JT ground state favors electron spin pairing, gener- 
ally leading to larger energy gains for low-spin states than 
for high-spin states. Molecular Hund-rule exchange and 
JT are therefore competing effects from the point of view 



of magnetism. Past studies indicated that in Gg(7 ions the 
JT interaction is in fact somewhat stronger than Goulomb 
exchange, leading to low-spin ground states OMZIIHIEI • ^ 
Nonetheless, the balance between JT and Hund's rule is 
to some degree compound dependent, leading to a close 
match in certain chemical environments. Spin gaps of the 
order of 100 meV between the S = ground state and the 
S = 1 excitation of the C^q or equivalently of the CqqIOu 
were reported in NMR studies of compounds like K4G60 
and Na2G6o, but a narrower 10 meV gap has been sug- 
gested for the n = — 2 fluctuating charge state in GsGeo 
[H] . The magnetism of TDAE-Geo pitTl) can also be at- 



^ Low-spin ground state is also observed in a number of di- 
radicals |10llll . but the situation there is quite different from 
Ceo ions. In ideal icosahedral Ceo ions the degenerate orbitals 
require a regular first Hund rule with maximum spin. Density- 
functional theory (DFT) calculations do of course confirm that. 
A low-spin ground state may be eventually obtained if, upon 
lowering the icosahedral symmetry by means of a JT distortion, 
the corresponding energy gain happened to be large enough to 
reverse the undistorted high-spin situation. In the aromatic di- 
radicals instead, two distant unpaired spins interact somewhat 
weakly through the molecular backbone, and a Heitler-London 
singlet ground state appears to be achieved, no distortions in- 
volved. 
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tributed to a spin triplet in charge fluctuating CgQ states 

m 

The main question which we wish to address here is 
what ground-state spin is to be expected for positive Cg,^ 
ions. The calculations will be done in parallel for Cg,^ and 
Cqq so as to emphasize the analogies and the differences 
that will emerge. Zero-point effects are also approximately 
included, and their effect is shown to be non negligible. 

Data on CgQ*" molecular ions do not appear to be read- 
ily available. In particular the published photoemission 
spectra of Cgo [TBItT7| imply Cgg final states. This n = 1 
state is affected by JT but clearly not by Hund's rule. 
Some charge-transfer compounds contain nominal Cqq ions 
While it is presently unclear whether higher charge states 
are accessible, acceptor compounds with nominal (AsF6)2 
Cgo and (SbFg)2C6o stoichiometrics have been described, 
containing again Cg^j" molecular ions. Recent data, includ- 
ing anomalously short spin-lattice relaxation times .18 , 
strongly suggest the possibihty of magnetism of the CgQ 
ions. 

Here we address the question of magnetism of Cgg*" ions 
by addressing quantitatively the competition between JT 
and Coulomb exchange in these ions. Both JT and intra- 
molecular exchange terms are calculated to be individually 
stronger in CgQ*" than in Cgg" (5yl9) . Preliminary calcula- 
tions and estimates based on the adiabatic approximation 
|S] did foreshadow the predominance of Hund's rule ex- 
change in Cgo'", as opposed to the predominance of JT 
in Cqq. In the adiabatic approximation however, the JT 
effect is postulated to be static, neglecting zero-point mo- 
tion of the carbon nuclei (responsible for turning JT from 
static to dynamic, even at T = 0). This simplification 
might be unsafe, as it is known to underestimate the JT 
energy gain ^8^. In turn, it might seriously impair our un- 
derstanding of the competition between high spin and low 
spin, through the neglect of very large corrections asso- 
ciated to the different zero-point vibrational quantum ki- 
netic energy associated with states of different spin. Here 
we estimate these zero-point effects and include them in 
the coupled JT-Hund's rule problem, thus correcting the 
adiabatic result. An accurate evaluation of this correction 
is in fact quite generally a formidable task. We show here 
that one lucky feature of Cgg*" ions is that their strong e- 
V coupling makes the basic zero-point correction a good 
approximation to the true result. Our bottom- line con- 
clusion will be that whereas quantum corrections to the 
ground state energy are confirmed to be large, they do 
not reverse the ground-state spin state for any of the Gqq 
ions. In particular the prediction that Cqq should be mag- 
netic is maintained. It is in fact reinforced, with the S* = 1 
ground state about 30 meV lower than the lowest S = 
state. 

This prediction does not of course imply that com- 
pounds containing C^q will by necessity be magnetic. 
Electron kinetic energy associated with hopping between 
molecules do in principle favor band electron spin pair- 
ing, and a standard nonmagnetic metallic state. However, 
any insulating states that could arise due to very nar- 
row bands, such as a Mott or more probably a Mott-JT 



insulator |2()|. is very likely to be magnetic. In this regard 
it seems interesting that both (AsFg)2C6o and (SbFg)2C6o 
acceptor intercalated materials, nominally containing dou- 
bly positive fuUerene ions were indeed found to be elec- 
trical insulators, with activation gaps of 0.22 and 0.64 eV 
respectively PT] . 

The starting point of this paper will be the static, adi- 
abatic JT-distorted state of all Cg^^ ions treated in the 
simplest model Hamiltonians that includes both e-v and 
e-e interactions. Each of these static JT wells, or "val- 
leys" , is characterized by a reduction of symmetry from 
icosahedral to some subgroup. New vibrational frequencies 
arise at each such valley. We determine these frequencies 
by evaluation of the Hessian energy matrix at the energy 
minimum 22 . We use the lowering of the zero-point vibra- 
tional energies from the undistortcd to the JT distorted 
state to estimate the leading quantum correction to the 
adiabatic approximation. The next-order quantum correc- 
tion would arise from weak tunneling between equivalent 
valleys. This is the conceptually important step that leads 
from static to dynamic JT, with full restoration of the 
undistorted icosahedral symmetry. The associated energy 
correction expected is however relatively minor in Gqq, 
in view of the large e-v couplings. In the limit of infinite 
coupling the tunneling corrections vanish, and all quan- 
tum effects coincide exactly with the zero-point lowering. 
Therefore, a direct quantitative check that the tunneling 
corrections are reasonably small already in the weakest 
coupling case of CgQ , and already negligible in Cgg justi- 
fies us in neglecting them for all other charge states. The 
accuracy of this neglect is, it should be noted, particularly 
good for the positive ions CgQ*", where the e-v coupling is 
larger than in Cgg" . 

This paper is organized as follows. Sect. El introduces 
the model and the parameters used in this calculation, 
which is then described in Sect. 13 along with the proper- 
ties of the JT valleys for all values n and S. The zero-point 
non-adiabatic corrections are described in Sect.^ and the 
overall results are finally discussed in Sect.|3 



2 The model Hamiltonian 

We begin by reviewing here the model Hamiltonian pre- 
viously introduced in Ref. p] to describe the electron- 
vibration coupling and Coulomb exchange of holes in the 
hu fivefold-degenerate highest occupied molecular orbital 
(HOMO), and of electrons in the threefold-degenerate ii„ 
lowest unoccupied molecular orbital (LUMO) of Cgo: 

H = Ho + Tvib + Kib + Hc-v + H^-c (1) 

where 

Ho ^ <^^cl^Cam, (2) 

am 
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Kib = E -^Q'a,, (4) 

iAfi 

-'-'e-v — / , 2 / ^ ^mm' ^^^A* '^crm^o- -m'AtJ; 

r iyl (7mm' f£ 

Hc-c = X! U;<^,£7'("^,m';fc,fc')4m4'm'C<T'fc'C<TA;(-6) 

are respectively the one-electron Hamiltonian, the vibron 
kinetic energy, the harmonic restoring potential toward 
the equilibrium configuration of neutral Cgo, the electron- 
vibron coupling (in the linear JT approximation) '19'"23J, 
and finally the mutual Coulomb matrix element represent- 
ing intra- molecular repulsion between the holes/electrons 

Here denote creation operators of either a hole 
in the HOMO or an electron in the LUMO, described by 
the single-particle wave functions (pmo-Cr). cr indicates the 
spin projection; m labels the component within the degen- 
erate electronic multiplets, according to the C5 character 
from the Ih D D5 D C5 group chain |2;-{ll24j . i enumer- 
ates the vibration modes of symmetry A {2 Ag, 8 Hh, 
and 6 Gg modes, the latter being JT active in the hole 
case only). CJJ^, are Clebsch-Gordan coefficients |22 of 
the icosahedral group Ih, for coupling hu/tiu states to 
phonons of symmetry vl. r is a multiplicity label, relevant 
for hu holes and Hg vibrations only, where it takes two 
values, 1 and 2 [HEII- QiA^i are the dimensionless molec- 
ular normal-mode vibration coordinates (measured from 
the adiabatic equilibrium configuration of neutral Ceo, in 
units of the length scale xo{uJi\) = ^hj (uji\ mc) associ- 
ated with each harmonic oscillator where mc is the mass 
of the C atom), and PiAfi the corresponding conjugate mo- 
menta. Finally, spin-orbit, exceedingly small in Ceo |25| . 
is neglected throughout. 

For holes, the JT model defined by Eq. ©, Q, ©, 
for Cqq is conventionally denoted a,s h" ® {A + G + H) : 
/i" refers to the hole occupancy of the hu HOMO, and 
A, G, H refer to the 2 nondegenerate Ag, 6 fourfold- 
degenerate Gg and 8 fivefold-degenerate Hg molecular vi- 
bration modes that are linearly coupled to in icosahe- 
dral symmetry |19II23II26| . For electrons, the JT model is 
t" ^ {A + H), referring to n electrons occupying the 
tiu LUMO, linearly coupled to the 2 Ag and 8 Hg vi- 
brational modes only. In all calculations we shall adopt 
the numerical values of the e-v coupling parameters gl^, 
listed in Tabled previously obtained from first-principles 
Density DFT electronic structure calculations in Ref. 19 . 
A recent DFT calculation [S^ based on a different func- 
tional reported couplings that are similar on the whole to 
those of Tabled the main difference concerning a closer 
competition between D^h and D^d valleys for the distor- 
tion of Cgo ■ The JT stabilization energy based on the 
DFT parameters of Ref. ^5] is only about one fifth of 
that found based on the intermediate neglect of differen- 
tial overlap (INDO) model |2Hj. We think that these ear- 
lier calculations, as well as more recent Hartree-Fock (HF) 
estimates [21] also suggesting large e-v couplings and en- 
ergy gains are somewhat less dependable. We believe the 



Table 1. Computed vibrational eigenfrequencies, and e-v lin- 
ear coupling parameters for the hu HOMO and tiu LUMO in 
Ceo Ei- 
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Hg 

261 


32.4 


3.0417 


-0.0045 


0.4117 


429 


53.2 


1.0587 


0.6131 


0.4886 


718 


89.0 


0.0103 


0.9950 


0.3500 


785 


97.3 


0.7836 


-0.0309 


0.2238 


1119 


138.7 


0.0514 


0.2151 


0.1930 


1275 


158.0 


0.4586 


0.2440 


0.1382 


1456 


180.5 


0.8482 


0.4530 


0.3152 


1588 


196.9 


0.7436 


-0.4488 


0.2893 



Table 2. The Coulomb parameters for Cg,^, defining ^c-o 
through JSJ and ||7|l, as obtained from the DFT calculations of 
Ref. (Sj. 





Parameter 


Value 
[meV] 


HOMO (C^o^) 








Fi 


15646 




F2 


105 




Fa 


155 




Fi 


47 




Fs 







U 


3097 


LUMO (c;7) 








J 


32 




U 


3069 



smaller DFT JT gains more realistic for two main reasons: 
(i) the HF calculations miss an important loss of correla- 
tion energy due to opening of the JT gap; (ii) preliminary 
results indicate that the HOMO photoemission spectrum 
based on the DFT parameters [21] is in very good agree- 
ment with experiment |16II17| . By contrast, for C^q, DFT 
computed electron-phonon couplings |!illl82IIH8ll84j appear 
to be significantly smaller than those obtained from pho- 
toemission of C^g |35| . While we can offer no explanation 
for that discrepancy in negative ions, we note that the 
DFT-derived couplings used here for C^^ are probably 
also slightly underestimated, however probably only by 
some 10 to 20%. 

The e-v couplings g^^ in Eq. |(SJ) are dimensionless, 
measured in the units of the corresponding harmonic vi- 
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brational energy quantum hiViA- Modes (e.g. the second 
Gg mode) with g^^ <^ 1 are thus only weakly coupled; 
conversely, modes (e.g. the lowest Hg mode) with a large 
gly^ > 1 are strongly e-v coupled. The numerical factors 

= 5*, k^<' = (5/4) ^ fc^f = 1 for the HOMO, and 
k^s = 32^ k^9 = 52 for the LUMO are introduced for 
compatibility with the normalization of the e-v parame- 
ters in Ref. . 

The Coulomb matrix elements are defined by: 



of that molecular ion Cq^, such that the average energy 
for each n 

i;-™(n) -en + C/^^^^^^-^. (11) 

We observe that U differs from a more common defini- 
tion of the Hubbard repulsion in lattice models, which 
involves the lowest state in each n-configuration: J7" 



E'^"^{n + 1) + E" 



w„,^-(m,m';fc,fc') =ydVydV (7) 

V*ma (r) i^') (r, r') Lfka (r) ifk'a' W) 

where u^.a' (r, r') is an effective Coulomb repulsion, screened 
by all other electrons of the molecule (eventually also 
by electrons in all other molecules in a solid state com- 
pound; but we shall focus here on the isolated ion). A 
detailed symmetry analysis ,5 shows that, assuming spin- 
independence of the orbitals, this set of coefhcients can be 
expressed as 



W(t,(t' {m, m'; k, k') 



r.r' ,A 



(8) 



in terms of a minimal set of independent Slater-type pa- 
rameters F"^'^ 1321 ■ A DFT estimate of these parameters 
was previously obtained in Ref. 5 , and for our calculation 
we adopt those values of the Coulomb parameters. They 
are reproduced for completeness in Table |21 Other sets 
of couplings, obtained by means of a simple (and clever) 
model for the HOMO and LUMO orbitals |S7| and by 
a fit to multi-configuration HF calculations [221 1 are in 
substantial agreement with each other, but they are both 
much larger than the DFT couplings of Ref. 0, due to 
complete neglect of screening. We believe that the actual 
Coulomb parameters of Cgo he somewhere in between the 
DFT couplings used here and the "bare" ones of Ref. |37l 
I29| . but most likely closer to the DFT ones, due to the 
large polarizability of Cgo. 

For the HOMO Coulomb parameters we use the short- 
hands 



Fi = F^o , F2 = F'^" 



F.^F 



1,1, -ffo 



(9) 



Fa = F 



2,2. Ha 



F.^F 



The Coulomb energy connected with the total molecular 
charge fluctuation in the HOMO, conventionally called the 
hole "Hubbard U" is given by the combination 



U = 



Fi 
5 



4^2 

45 



El 
9 



F\ 
9 



(10) 



For electrons in the LUMO, we have instead an electron 
Hubbard U given by C/ = i^^f /3-F^f /3 and a Hund-rule 
exchange J — F^^ /2. 

In either case, of electron or of holes, U defines an 
average Coulomb repulsion within the multiplet of states 



1) - 2£;'"™(n). This second defini- 



tion is inconvenient here, since it depends wildly on n. The 
Hubbard U's are given here for completeness, and because 
they will be useful in different contexts. However we must 
recall that the U term is irrelevant for the determination 



of the ground-state spin of each Cg^*^ ion 



3 The adiabatic JT valleys 



Though a rather idealized description of real CgQ ions, 
the model Hamiltonian (^), does nevertheless not lend 
itself to an exact solution. The vibronic eigenstates are 
generally complicated combinations in the direct product 
of each (up to 252-dimensional) fixed-n electronic space, 
times the infinite-dimensional space of the vibrational de- 
grees of freedom. Even with the help of spin and orbital 
symmetries, exact solutions of the quantum problem 
are only available in two limiting cases, namely the limit of 
weak e-v coupling '5"'B 5II5^ and that of infinitely strong 
e-v coupling [7,22,26,401. For Cgi^^, where the couplings 
range from intermediate to large, some approximations are 
therefore called for. If the phonon kinetic term Tvib is ne- 
glected in the so-called adiabatic approximation, the dis- 
tortion operators QiAii are replaced by c-number coordi- 
nates. This approximation yields the leading ground-state 
energy lowering c>c '^{glAf'^^'^iA, exact in the limit of large 
e-v couplings 00. [Note however that, in this limit, 

both initial assumptions of harmonic vibrations I^J and 
of linear e-v coupling jnj become anyway questionable]. 
In Sect. 21 we shall deal with the leading quantum correc- 
tions to the adiabatic approximation oc J2i9iA)^^'^ij^ t)y 
taking zero-point energy shifts into account. 

The classical treatment of the vibration coordinates 
breaks the full molecular symmetry (here icosahedral sym- 
metry) in all configurations possessing a nonzero distor- 
tion QiAfi {A = Gg,Hg). Therefore states of different 
icosahedral symmetry representations are in general inter- 
mixed, leaving only the total number of holes n, the total 
spin S and its projection Sz conserved in the adiabatic 
ground state. Here we will assume the orbitals to remain 
basically unchanged upon JT distortion, thus neglecting 
any small JT-induced change of the Coulomb Hamiltonian 
Hc-c- The latter is therefore still determined according to 
Eqs. H6I7I8I) by the same parameters Fi of Table [21 as in 
the undistorted icosahedral configuration. We will also as- 
sume no change of the vibration frequencies and cou- 
plings upon charging. While this is of course at vari- 
ance with much established evidence showing vibration 
frequency shifts in the percent range per each added elec- 
tron, it is perfectly in line with the other approximations 
intrinsic in our model, and in its solution. We leave the 
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Ag modes out of the adiabatic calculation, since despite 
their nonzero linear e-v coupling they simply contribute a 
trivial spin- and symmetry-independent term 

i 

a = 1.79 meV (HOMO), a = 2.90 meV (LUMO), 

to the total energy E^^°{n) [Eq. ^(TT^ ]. which could effec- 
tively be included into U. Because of particle-hole sym- 
metry (exchanging creation and annihilation operators) 
of the Hamiltonian (excluding Hq , the Ag modes and the 
average U contribution), positive charges n > 5 can al- 
ways be mapped onto n < 5, and negative n > 3 onto 
n < 3. 

For each charge n and spin S and as a function of 
the vibration coordinates Q, the lowest adiabatic poten- 
tial surface ^^^(Q) is obtained as the sum of the lowest 

eigenvalue of 7Je_e+J?e-v(Q) in the n-electron spin-S* sec- 
tor, plus the harmonic restoring term Vvib(Q)- Allowing 
the 64 (6 X 4 Gg ^ plus 8x5 Hg.) vibration coordinates Q 
to relax, we determine the optimal distortions Qmin by full 
minimization of l^s(Q) in the space of all vibration co- 
ordinates Q. Details of this minimization were previously 
reported in Ref. [KIHT] . 

Table |31 summarizes the optimal adiabatic energies 
Eadicihin, S) — t^^5(Qmin)- The Separate vibration, e-v, 
and exchange contributions are specified. The results show 
that the adiabatic valley minima of positive Cgo ions are 
systematically lower for high-spin states, while the adia- 
batic outcome for negative ions favors marginally low spin, 
and is essentially uncertain, with small and poorly reliable 
energy differences between high- and low-spin states. 

Positive and negative ions differ qualitatively also in 
the geometric nature of the JT distortions characterizing 
the minimal valleys in the space of coordinates Q. The 
H model applicable to electrons in CgQ" leads to a 
flat continuous manifold of equivalent points. The valleys 
are really troughs, which could be imagined as "Mexican 
hats" [7ll42n43| : the study of these 2-dimensional (n = 1, 2) 
and 3-dimensional (n = 3) troughs is a classic topic in JT 
physics PIUI] . 

Conversely, for holes in C6o+ the h" ® {G + H) JT 
model gives rise to discrete sets of isolated minima [T^ll^ 
I26ll45j . The nature and symmetry of these discrete valleys 
depends on the details of the e-v couplings and on the 
interplay with Coulomb interaction. We shall describe here 
in some detail the minima obtained in this calculation. 

To identify these minima, we generate about a hundred 
randomly distributed distortions away from the Ih high- 
symmetry point, and let the vibration coordinates relax 
from there to the closest energy minimum, by combined 
standard (simplex and conjugate-gradients) minimization 
algorithms. We then apply the symmetry operations of 
the icosahedral group to the each minimum found, locat- 
ing all possible equivalent minima. Although the method 

^ The Gg modes are irrelevant for the negative ions, but rel- 
evant for positive ions 



employed is not exhaustive, the application of symmetry 
and the thoroughness of our search, carried out with a 
large variety of starting points guarantees in practice the 
retrieval of all relevant minima. We can readily discard, 
for example, the few cases where minimization led to sad- 
dle points or to secondary, local minima, based on simple 
comparison of the energy values. In this way, for each n 
and 5*, we finally obtain the set of equivalent valleys that 
are global adiabatic energy minima for Cqq*" . The amount 
of "radial" distortion at the adiabatic minima for each 
mode is tabulated in the Appendix. 

In Table ^ we summarize some global properties of 
the JT energy valleys for Cg^*^ in all spin sectors. In these 
many-mode JT systems, the local symmetry of an optimal 
distortion is described in terms of the subgroup Giocai C Ih 
of symmetry operations which leave that minimum invari- 
ant. The minima in the simple n — 1 S — ^ case, where e-e 
interactions are irrelevant, are believed to be six valleys of 
D^d symmetry |T^1I^S1I^ - the alternative possibility of 
ten D^d valleys apparently disfavored by the specific cou- 
plings obtained for fuUerene 19 . Table ^ reveals a sym- 
metry between the minima for n — A S — 2 and n = 1 
5 = ^, as well as for n = 2 S" = 1 and n = 3 5' = |. 
These symmetries are at first sight surprising, but can be 
readily explained by applying a particle-hole transforma- 
tion to the fermion operators of only one spin kind, in the 
fully spin-polarized states. This transformation maps the 
Hamiltonian matrices of an n-particle states (apart from 
a constant exchange term) into those of (5 — n) particle 
states, with a sign change of the vibron interaction -ffo-v 
that shifts each minimum Qminto the opposite locations 
-Qmin- We have therefore y4^^(Q) = V^^'^^^i-Q) + C, 

with C = -f-Fa - |i^3 - |F4 = -238 meY. The same 
connection between the n = 2 S = 1 and the n = 3 
5' = I adiabatic potentials involves an energy shift C = 
— |-p2 — ^-^3 — — nieV. The same symmetry re- 
lates n — 5S— ^ton = OS = 0, where no JT distortion 
takes place. 

In general, the number of JT valleys (third column of 
Table 2J| is related to the local symmetry. It is basically 
given by the ratio |/h|/|Giocai| of the order of the icosa- 
hedral group (120) to the order of the invariant subgroup 
Giocai- However, n — 5 holes make an exception to this 
rule. Here, at half filling, in addition to the Ih symme- 
try, the system is particle-hole symmetric, i.e. invariant 
under exchange of creation and annihilation operators of 
both spin kinds. This transformation leaves the Coulomb 
Hamiltonian -ffo-c invariant, and again changes sign of 
the vibronic interaction Hc-v'. hence, given any minimum 
Qmin, its opposite —Qmin is also an equivalent minimum 
of the potential energy surface. For S = ^, this leads to 
a doubling of the minima: the local C2h symmetry would 
lead to 30 minima, but 30 additional equivalent minima 
are added in the opposite positions by particle-hole sym- 
metry. For S — ^ instead, the number of minima remains 
30, since for each minimum there is one of the Ih symmetry 
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Table 3. The three contributions (restoring harmonic potential, e-v couphng, and e-e exchange - the [Un{n — l)/2] term is 
excluded from He-e) in the adiabatic energy iJadiab for Cg,^ in all possible charge and spin states. As usually occurs in JT, 
the contribution ^//o-v) = —2(Vvib) = — 2Vvib(Qmin)- The following columns report the leading zero-point energy correction 
Szcro [Eq. I|15|l ] and the sum of this correction to the adiabatic energy, representing our best estimate of the energy of each 
valley, relative to the undistorted state of the ion. For C^q, the last column reports the exact energy obtained by Lanczos 
diagonalization. All energies are in meV. Boldface indicates the ground-state total energy at each n (high spin for cations, low 
spin for anions). For positive ions, comparison with the exact result in the last column indicates a fairly good approximation 
already for n =1; and we expect the accuracy to be even better for n = 2,3,4,5 holes. 



n S 








-Eadiab 




-£'adiab ~l~ -E'zcro 


-E'cxact 


1 1/2 


69 


-138 





69 


-51 


120 


-103 


2 


270 


-540 


141 


-129 


-72 


-201 




1 


99 


-197 


-43 


-142 


-92 


-234 




3 1/2 


267 


—534 


99 


— 168 


— 122 


—290 




3/2 


99 


-197 


-123 


222 


-92 


-314 




4 


361 


-723 


162 


-200 


-133 


-332 




1 


229 


-459 


19 


-211 


-134 


-345 




2 


69 


-138 


-238 


308 


-51 


-359 




5 1/2 


308 


-616 


105 


-203 


-150 


-353 




3/2 


169 


-338 


-87 


-256 


-98 


-354 




5/2 








-397 


397 





397 




1 1/2 


38 


-77 





38 


-113 


-152 


-76 


2 


149 


-298 


56 


93 


-115 


-207 




1 


38 


-77 


-32 


-71 


-113 


-184 




3 1/2 


113 


-225 


28 


-85 


-171 


-256 




3/2 








-97 


-97 





-97 





operations, a C2 rotation, that transforms this minimum 
into its opposite point'^. 

Table 01 contains also some information about the con- 
nectivity of the minima in Q space. In many cases, the 
specification of the number of first, second, etc. neighbors 
of a given minimum is sufficient to clarify completely the 
topology of the minima in the 64-dimensional space. In 
particular, the wells of the n = 1 S* = i, of the n — 2 
5 = 0, and of the n = 4 S = 2 surfaces are located on 
the six vertexes of five-dimensional regular simplexes, gen- 
eralizations of the 3D tetrahedron, each minimum being 
equidistant to all the others. In analogy, the connectivity 
of the 10 Z?3d minima for n = 4 S' = is the same as 
that depicted in Fig. lb of Ref. [221 for a different situa- 
tion. For the other cases of lower symmetry, the number 
of neighbors of any given order must be complemented by 
some extra connectivity information, for which we refer to 
previous work j^H- We only observe that for 71 = 2 = 1 
(equivalently for n = 3 5 = |), each of the 15 minima is 
linked to four nearest-neighbor minima, which, in turn, are 



^ The real-space inversion, acting on the real-space displace- 
ments, leaves all JT-active normal-mode coordinates invariant, 
since they are even {g) under inversion. 



linked to more minima, forming a completely connected 
regular polytope. For n — 3 S = and n = 4 S" = 1, the 
30 minima are divided into 6 pentagonal "clusters" of five 
nearest- neighboring minima. In contrast, for n = 5 S ~ ^, 
nearest- neighbor wells come in pairs. Finally, the 30 C2h 
minima for n = 5 S = |, show the largest connectivity, 
and sit at the vertexes of a highly symmetric polytope. 

Figuren]pictures the electronic state of the holes at the 
adiabatic valleys of Cqq. As typical of JT, distortions are 
such that relatively large gaps open between empty and 
filled levels. For n > 1 this single-determinant picture is 
only approximate: it represents the dominant uncorrelated 
configuration within the fully correlated exact electronic 
ground state at each well. 

The adiabatic valleys described above are the joint re- 
sult of e-v couplings and of e-e Coulomb interactions. One 
can understand better the specific role of Coulomb inter- 
action for the JT distortion by comparing the number and 
symmetry of the valleys described above to those obtained 
by the same procedure in a hypothetical noninteracting 
case, obtained by setting all Fj = 0. Doing that, we find 
basically the same general picture of valleys for n = 1 (ob- 
viously), but also ioT n = 2, n = 3, n — A S — 1, 2 and 
n = 5 S* = |, |. In all remaining cases the e-e correlations 
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Table 4. Number and local symmetry of the JT minimal valleys for each charge n and spin S in Cg(f ions. The 5th column 
gives, for each given valley, the number of equivalent valleys that are first, second, etc. neighbors of that valley in Q space. The 
last column gives the total magnitude of the dimensionless JT distortion at each minimum. 
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22144462 




2 
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1.58 


5 1/2 


60 
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1 2 


244422242 242224442211 


3.27 


3/2 


30 




C2h 


8 12 8 1 


2.46 


5/2 


1 




Ih 
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1 1/2 


oo 




or higher 


oo 




0.91 


2 


oo 




or higher 


oo 




1.79 


1 


oo 




or higher 


oo 




0.91 


3 1/2 


oo 




or higher 


oo 




1.55 


3/2 


1 




Ih 











also introduces a qualitative change in the topology and 
local symmetry of the minima. Without interactions, 15 
D2h minima would replace the 10 D^d minima for n — A 
S* = 0; and 20 D^d minima replace the 60 C2h minima 
for n = 5 5 = 1/2. In all cases (except of course n = 1 
and the symmetric n = 4 S = 2), the amount of distor- 
tion is a few per cent larger in the uncorrclated case than 
when exchange is included, as expected from competing 
interactions. For example, without correlation, the total 
distortion for n = 2 5 — would be 3.16, i.e. exactly 
twice the distortion for n = 1, as routinely occurs in JT 
when exchange plays no role. 



4 Zero-point quantum corrections 

The vibrational quantum kinetic energy, a quantity ne- 
glected at the adiabatic level, does in fact contribute quite 
substantially to the JT energetics of these multi-mode sys- 
tems characterized by intermediate couplings. The zero- 
point energy gain associated to the softening of the vibra- 
tional frequency at the JT-distorted minima is illustrated 
in Fig. 121 It represents the leading quantum correction to 
the static JT energetics 22,41,46 . To compute this cor- 
rection for each charge n and spin 5, by finite differences 
we evaluate the Hessian matrix of the second-order deriva- 
tives of the lowest adiabatic potential sheet, at each of the 




^adiab+'^zero 



distortion 



Fig. 2. A pictorial illustrating the origin of the zero-point 
energy gain due to vibrational softening at the JT well. The 
horizontal axis represents a generic distortion coordinate, the 
horizontal thick lines represent the quantum ground states in 
the non-JT and JT potential wells. 



adiabatic JT minima: 



H{n, S)aa' 



dQa dQa' 



(13) 



Qmi„(",S) 



where a is a collective index for {i A The normal-mode 
frequencies ujj in each well are computed by taking the 
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E [meV] 



-100 



-200 



n=l 



n=2 



n=3 



S= 1/2 



S= 



S= 1 



S= 1/2 



S= 3/2 



E [meV] 



-100 ■ 



-200 ■ 



n=4 



s= 



S= 1 



S= 2 



S= 1/2 



n=5 

S= 3/2 



Fig. 1. A line-spectrum representation of the dominant single-determinant electronic state in the optimal adiabatic configuration 
of CgQ*" for each spin 5*. Levels are ordered for holes. 



square roots of the eigenvalues of the dynamical matrix 
defined by 



S)aa' 



- w, 



1/2 



[U{n,S)a, 



>UJa\ iO, 



1/2 
a' ' 

(14) 

The —Saa'^a term removes the restoring term from V^jb in 
V^%, and the LoU"^ are introduced by the standard change 
of variables to correct for the different "mass" coefficients 
of different coordinates Qa in Tvib- By retaining the har- 
monic expansion of the adiabatic potential around a min- 
imum the quantum ground state energy of this potential 
well is estimated at ^fiCjj above the classical minimum, 
due to zero-point motion. As illustrated in Fig. |21 the dif- 
ference between this and the original zero-point energy 
'Ylia \^^a at the neutral-molecule harmonic minimum pro- 
vides the leading quantum correction 



j iAfj. 



iA 



(15) 



to the "classical" valley energy i^adiab- Ezcio{n, S) in JT 
problems is systematically found to be negative in sign, 
corresponding to softer vibrations, and a shallower valley 
bottom at the new minima than for the n = non-JT 



molecule. Table 01 reports the zero-point correction com- 
puted as described above, and the estimate of the ground- 
state energy obtained by adding the correction Ezcm to 

-^adiab • 

For Cgo" , the zero-point corrections favor low-spin states 
for both n = 2 (already a, S = ground state at the adi- 
abatic level) and n = 3, where the zero-point correction 
only affects S ~ 1/2, lowering it well below the uncor- 
rected S = 3/2. Note however that in the t" « 77 JT for 
the anions, the two (n = 1 and 2) and three (n — 3) 
"soft" modes along the trough are associated to vanish- 
ing frequencies ujj. They formally contribute a vanishing 
zero-point energy to the first sum of Eq. H15|l . In our cal- 
culation this produces unreasonably large (negative) zero- 
point corrections, larger even than the adiabatic energies. 
These null terms are correct only in the limit of infinite 
coupling, where the "size" of the fiat trough is infinite, and 
indeed the free pseudorotational motion carries no zero- 
point energy. For Cg,^, where coupling is finite, and not 
especially large, the trough has a finite size of order g'^, 
which provides some amount of quantum confinement, as- 
sociated to a zero-point kinetic energy of order as was 
described in Ref. [7li43) . These extra "confinement" cor- 
rections should be especially sizable when the vibrational 
wavefunctions have the nontrivial nodal structure due to 
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the boundary conditions associated to an electronic Berry 
phase, i.e. for n = 1, n = 2 S* = 1, and n = i S — 1/2 
but should not change the conclusion that all Cg,^ ions 
should favor low spin. 

To get an estimate of the validity of the approxima- 
tions employed here, we can compare with the ground 
state energy of Cgg obtained by means of an essentially ex- 
act calculation (a Lanczos diagonalization on a truncated, 
but well converged, basis) of our Hamiltonian ^ with 
the same parameters. The exact energy gain obtained is 
—76 meV, which falls in between the adiabatic (—38 meV) 
and zero-point corrected (—152 meV) values. Not surpris- 
ingly, in this rather weakly coupled case, the zero-point 
correction largely overshoots the correct value. As charges 
n > 1 are associated to larger distortions, thus effectively 
to stronger couplings, we expect that the zero-point cor- 
rected energy should be of better quality there. However, 
the relative size of the zero-point and confinement cor- 
rections to the adiabatic energy makes the whole strong- 
coupling expansion rather questionable for all the anions. 
Thus the G^^ results as reported here serve mostly for 
comparison with those of the cations. For cations we ex- 
pect in fact the quality of our zero-point corrected results 
to be substantially better. 

In the Cg(^ ions too, in fact, the zero-point energy cor- 
rections are rather large, but here the magnitude of the 
adiabatic energy gains is even larger. The difference with 
the anions is due both to the larger couplings and to the 
structure of localized minima as opposed to a flat trough. 
Comparison of Table 13 and Table 0] shows that situations 
with a larger number of minima generally yield larger val- 
ues of the correction |i?zGio|- Quantum kinetic energy in 
this respect behaves a bit like statistical entropy, in that 
it favors numerous shallower minima against few deeper 
ones. The zero-point correction peaks at « —150 meV for 
n — b S = ^, where the minima are shallow and the low- 
est vibrational frequency is as small as Cji ~ 10 meV in 
this case. On the contrary, few well spaced minima, as 
for rt = 1 5 = 1/2 and n = 2 5 = are associated to 
a smaller zero-point gain of the order « —100 meV. In 
the close competition between Coulomb physics (Hund's 
rules) and JT physics (anti-Hund behavior), the zero-point 
correction is not irrelevant. As shown by the last column 
of Table 13 it reduces drastically the large adiabatic spin 
gap between the high-spin ground state and the lowest 
spin excitations for n = 4 and for n — 5. Remarkably, the 
zero-point correction favors the high-spin state instead in 
Cqq: the excitation energy to S* = is enhanced by the 
quantum correction from 13 meV to 33 meV. Thus al- 
though the estimated spin gap is not especially large nor 
especially reliable, this circumstance makes in our view 
the prediction of a S* = 1 magnetic C^q ion stronger than 
for the other cations. 

Could we go beyond the zero point correction, and get 
more accurate results? Not easily at this stage. The zero- 
point correction represents the g'^ term of a large-coupling 
expansion, where the adiabatic energy £^adiab is the lead- 
ing (g^) term. The next corrections to be considered, of 
order g~^ and higher, are associated with anharmonicity 



of the valleys and tunneling amongst them (in the holes 
case), and also to mixing of the upper adiabatic poten- 
tial surfaces with the associated geometric-phase effects 
[7ll23ll47ll48| . A full quantitative description of these ef- 
fects would imply a much more sophisticate treatment of 
the quantum problem than the simple semiclassic expan- 
sion applied here, and that is beyond the scope of the 
present paper. To get an estimate of the importance of 
these higher-order corrections, we can again compare, now 
for Gqq the adiabatic (—69 meV) and the zero-point cor- 
rected (—120 meV) ground-state energy gains to the exact 
on, obtained by a well converged Lanczos diagonalization, 
which is —103 meV. We see that although E^cm still over- 
shoots the quantum corrections the adiabatic gain, the 
residual error of 17 meV is now fairly small. The accu- 
racy expected for the zero-point corrected energy gains 
for all other cations Cgj^^^"*" is even better. In fact, the 
distortions there are larger, tunneling is suppressed, and 
the approximation of isolated harmonic JT valleys should 
be much better than for n = 1. Thus, even if a check with 
the Lanczos method would be too cumbersome here be- 
cause of the excessive basis size required for ti > 1, our 
zero-point corrected results of Table 01 should be regarded 
as practically exact, for the assumed coupling parameters. 

5 Discussion and Conclusions 

The main output of this paper is a determination of the 
ground state spin, energy, and distortion magnitude of JT 
distorted Cg,^ molecular ions. 

Calculations include both e-e and e-v interactions as 
derived from earlier first-principles calculations, as well as 
very important quantum vibrational effects, due to the 
small mass of the carbon nuclei. The latter are taken into 
account approximately by including the changes of vibra- 
tional zero-point energy from the undistorted to the dif- 
ferent distorted states. These zero-point corrections are 
shown to be generally large. In Cg,^ they are even capable 
of turning a high-spin ground state into a low-spin one. 
As the coupling and thus the distortions are fairly large, 
quadratic and higher-order (in Q) e-v interactions and vi- 
bration anharmonicity could also be relevant. The present 
calculation was however carried out strictly in the linear 
e-v coupling approximation. 

The parameters used in this calculation, both for e- 
e and e-v interaction could be somewhat underestimated 
by the local density approximation used in their deter- 
mination, as discussed in Rcf. 5,19 . Consequently, both 
the Coulomb repulsion and the JT effective e-e attraction 
calculated within the local density approximation might 
need some correcting in their absolute values. The bal- 
ance between these two opposing interactions is delicate in 
Cg|7 ions (as demonstrated by the presence of both high- 
spin and low-spin local ground states in different chemical 
environments 9,15,49,50,51,52,53 ), but low-spin ground 
states are generally believed to prevail, in accord with the 
present calculation. Contrary to that, in Cgg*", Hund-rule 
magnetism and high-spin ground states are predicted to 
dominate. Experimentally, so far we found no evidence 
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concerning gas phase CgQ ions that we could usefully ad- 
dress. We do not know at this moment whether the lifetime 
of positive ions against fragmentation would permit exper- 
iments to be conducted, and if so whether the ground-state 
spins and JT distortions could be determined and com- 
pared with our predictions. Singly charged C^q ions have 
been created in solution j231j in molecular-beam photoe- 
mission |16II17| . and in storage rings (55j: but they are 
irrelevant to our present question. 

More encouragingly, higher positive charge states have 
been pursued in solid-state acceptor compounds (56k57l . 
The interplay of intra- molecular physics with electron hop- 
ping in a hypothetical fullerene-cation based solid-state 
conducting compound is very interesting. Our results sug- 
gest that magnetism should be important in these com- 
pounds, at least so long as the hole bands remain rela- 
tively narrow. The HOMO bands of the proposed accep- 
tor compounds (AsF6)2C6o and (SbFg)2C6o can indeed 
be expected to be rather narrow. In (AsF6)2C6o, based on 
a bet unit cell with lattice parameters a ^ b = 12.8 A, 
c = 12.4 A, we can extract a fuUerene-fuUerene inter- 
center distance 10.97 A, significantly larger than 10.02 A 
of pure Ceo. On the other hand, the on-site Hubbard U 
estimated for holes is comparable to or slightly larger than 
that for electrons. It seems reasonable therefore to surmise 
that hole-doped (AsF6)2C6o could be, like the electron- 
doped K3NH3C60 Mi, as well as A4C60 jSMil , a Mott-JT 
insulator. Unlike the electron-doped compounds, where 
the JT distorted Cg,^ molecular ion was in a low-spin con- 
figuration {S = 1/2 for n = 3, S' = for n = 2,4), we 
expect in the divalent acceptor fuUerene compound that 
the positive CgQ ions should be in a high-spin, S = 1 JT- 
distorted state. If so, they would constitute an exciting 
example of molecular Hund's rule magnetism. The mag- 
netic coupling between neighboring Cqq S — 1 ions in 
these compounds should most likely be weak and antifer- 
romagnetic . It might be interesting to pursue this and 
to seek a Neel state of some kind at very low temperatures. 
Pressure studies of the electron-doped compounds have re- 
vealed insulator-metal transitions such as that of Rb4C6o 
j(i2| . and also insulator-superconductor transitions such as 
that of K3NH3C60 Pressure studies of fuUerene ac- 
ceptors could be an interesting line to pursue in the future. 
More generally, the realization of a well characterized ful- 
lerene hole conductor remains a worthy challenge for the 
future. 
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Table 5. The dimensionless JT distortions \QiA\ at the min- 
ima of CgJ~ , for each Hg mode and allowed spin S. 



n 


S 


dimensionless distortions of Hg modes 


1 


1/2 


0.412 0.489 0.350 0.224 0.193 0.138 0.315 0.289 


2 





0.813 0.965 0.691 0.442 0.381 0.273 0.622 0.571 


2 


1 


0.412 0.489 0.350 0.224 0.193 0.138 0.315 0.289 


3 


1/2 


0.706 0.838 0.600 0.384 0.331 0.237 0.540 0.496 


3 


3/2 


0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 



A Appendix 

Tables El and report the individual amount of JT dis- 
tortion associated to each normal mode of Cgo at any 
equivalent adiabatic minimum of C^5^, for aU spin-S* sec- 
tors. The distortions are given in the dimensionless units 
of the coordinate operators The appropriate length 

scales xo(wiyi) — (lOiA nT-c) for the Gg and Hg modes 
of Ceo are: 76.3, 70.4, 60.3, 50.3, 46.1, 43.0, and 103.7, 
80.9, 62.6, 59.8, 50.1, 47.0, 43.9, 42.1 pm, respectively. 
For Cqq the amount of "radial" distortion is independent 
of the point chosen along the trough. As expected, the 
modes characterized by the strongest coupling show the 
largest distortion. The mode associated to the largest dis- 
tortion is therefore the lowest Hg mode (see Table ^1 in 
the cations. Note also that the Gg modes do not contribute 
to the distortions of -D5 symmetry, but contribute to all 
the lower-symmetry minima. 
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